Magnetic-field-induced superconductivity and superfluidity of W and Z bosons: 
in tandem transport and kaleidoscopic vortex states 
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We show that in a background of a sufficiently strong magnetic field the electroweak sector of 
the quantum vacuum exhibits superconducting and, simultaneously, superfiuid properties due to 
magnetic-field-induced condensation of, respectively, W and Z bosons. The phase transition to 
the "tandem" superconductor-superfiuid phase - which is weakly sensitive to the Higgs sector of 
the standard model - occurs at the critical magnetic field of 10^° Tesla. The superconductor- 
superfiuid phase has anisotropic transport properties as both charged and neutral superfiows may 
propagate only along the magnetic field axis. The ground state possesses an unusual "kaleidoscopic" 
structure made of a hexagonal lattice of superfiuid vortices superimposed on a triangular lattice of 
superconductor vortices. A weak electric field will induce both superconducting and, unexpectedly, 
superfiuid fiows. 

PACS numbers: 12.15.-y, 13.40.-f, 74.90.+n 
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It is know that an extremely high magnetic field of 
hadronic scale may lead to plenty of unusual effects both 
in (dense) matter and in the quantum vacuum. The chi- 
ral magnetic effect [1^ provides a particularly interesting 
example: a charge-parity-odd matter may generate an 
electric current along the axis of the magnetic field [3 . 
The corresponding conditions may be realized in noncen- 
tral heavy-ion collisions [2] in which a hot quark matter 
is created in a background of extremely high magnetic 
fields [H |5] . Similar conditions may have existed in the 
very early moments of our Universe [9 . The strong mag- 
netic field affects phases of the cold dense matter in the 
cores of strongly magnetized neutron stars [6 . 

Because of quantum effects an empty space may also 
exhibit quite unusual properties in a sufficiently strong 
magnetic background. In the background of a relatively 
low magnetic field of QED scale the vacuum should be- 
come become optically birefringent [7 . The hadron-scale 
magnetic field should lead to magnetic catalysis [8] , which 
implies, in particular, a steady enhancement of the chi- 
ral symmetry breaking in QCD vacuum as the external 
magnetic field strengthens. 

More recently it was found that the vacuum becomes 
an electromagnetic superconductor in sufficiently strong 
external magnetic fields JOl [IlJ • The superconductivity 
of, basically, empty space, is mediated via spontaneous 
creation of a (charged) p-meson condensate if the mag- 
netic field exceeds the critical value of 



^QCD ^ ^2/g _ ;^ Q ^ iQi^Xesla, 



(1) 



where Mp = 775.5 MeV is the mass of the p meson. The 
ground state of the vacuum superconductor is character- 
ized by inhomogeneous ground state of a very particular 
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geometric structure [12 , possessing intriguing metamate- 
rial ("perfect lens") properties [13]. The magnetic fields 
of the required strength ([l]) may be created on Earth in 
heavy- ion collisions at the Large Hadron Collider [5 . 

We show that as the background magnetic field 
strengthens further, the Standard Model experiences a 
second superconducting, and, simultaneously, superfiuid 
transition associated with a condensation of the W and 
Z bosons at a larger critical magnetic field: 



B, 



EW 



M^/e 



1.1 X 10^^ Tesla, 



(2) 



where Mw = 80.4 GeV is the mass of the W boson. 
The condensation of the W bosons at the magnetic 
field higher than the critical value ([2| was predicted by 
Ambj0rn and Olesen in Refs. [14 . The key idea here 
is that the vacuum of charged vector particles (i.e., of 
the W mesons) is unstable in the background of a suf- 
ficiently strong magnetic field provided these particles 
have anomalously large gyromagnetic ratio g^ = 2. The 
large value of g guarantees that the magnetic moment 
of such particles is too large to withstand a sponta- 
neous condensation at sufficiently strong external mag- 
netic fields. We show that the W condensate leads to 
exotic superconducting and superfiuid effects. 

The bosonic part of the electroweak sector of the Stan- 
dard Model is described by the Lagrangian, 
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where $ is the complex Higgs doublet which interacts 
with SU{2)l and u\l)x gauge fields (W^ and X^, re- 
spectively) via the covariant derivative D^ = d^ — 
igT^W^/2—ig'Xp/2^ and r^ are the Pauli matrices. The 
corresponding field strengths are IV^j, 
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FIG. 1. 

(c) the 
fields 



(a) The W condensate (IVKI^)"*^ , (b) the gain in the energy density ([6l due to the condensation, 6E = {E)w — {E)w^o, 
Higgs expectation value ($^$)^/^. Eq. ( |l7| , and (d) the superconductivity transport coefficient (20) vs. the magnetic 
at fixed values of the Higgs masses Mh — (1.1 .. . 2.0) M^ (all quantities are averaged over the transversal plane). 



The mexican hat potential in Eq. ([s]) breaks the elec- 
troweak symmetry down to the electromagnetic sub- 
group, SU{2)l X U{l)x -^ U{l)em because the Higgs 
field ^ acquires a quantum expectation value, {^) ^ 0. 
In the unitary gauge, {^) = {O^v)^^ the third com- 
ponent of the non-Abelian gauge field W^ mixes with 
the Abelian gauge field X^ providing us with the mas- 
sive Z^ boson and the massless electromagnetic field A^: 



W^ = smeA^ 



COS OZn^ X^i = cos OAl, 



sin OZ^ where 



is the electroweak mixing (Weinberg) angle with e = 
^sin^ = g' cosO being the electric charge. 
The classical equations of motion are: 

= d^W^^ + ge^'^'^W^^W^^-ig [(D^<l>) V'^^-Zi-c] /2 , 
= 5^ V - ig' {D^^)^^ - h.c) /2 , (4) 





-D^D^^ 



The instability of the vacuum in the presence of the suf- 
ficiently strong magnetic field was clearly demonstrated 
by Ambj0rn and Olesen in Ref. [14 , and we briefly re- 
peat their arguments here. Consider a uniform time-in- 
dependent magnetic field directed along the third axis, 
5ext,i = ^ext^i3 (wc always take eB > for simplicity of 
our notations). Then the quadratic part of the transverse 
(with respect to the magnetic field axis) components of 
the W^ = W- field in Eq. ^ reads 



sc^^\ = {wlM) 



ieBe^ct 



-ieBey 



Wi 
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(5) 



(with Mw = ^'^/2), while the mass terms of the lon- 
gitudinal components W3 and VFq, and of other vector 
particles are not affected at the classical level. The 
mass eigenvalues of Eq. (l5| are /i^ = M^ ± eB. One 
of the masses, /i_, vanishes at the critical value Be 
of the magnetic field ^. This mass becomes purely 
imaginary ai B > Bc^ thus signaling a tachyonic in- 
stability towards condensation of the transverse com- 
ponents of the W^ field, with the unstable eigenvector 
{Wi, W2) = {W, -iW)/2, where W is d. scalar field. 

Since we consider the solutions in the transverse 
(xi, X2) plane for the transverse components of the fields, 
it is very convenient to use the complex notations for 
the coordinates, z = xi -\- 1x2^ and for the vectors 



Oi + i02, O = Oi- i02 and O12 = -^{00 - dO). 
Notice that H/^t ^ w. 

We use a symmetric gauge for the external magnetic 
field, Aext,i = —BX2/2 and Aext,2 = Bxi/2^ so that the 
corresponding covariant derivative is Pext = d — ieAQ^i = 
9 + ez5ext/2. The W component of the W~ corresponds 
to the /i+ eigenvalue of the operator in Eq. ([5|, so that 
it is not condensed. Thus, we put I^ = as it does not 
lower the energy, and continue to work in the unitary 
gauge $ = (0, ^)^, where is a real- valued field. 

In complex notations the energy density is: 
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and the equations of motion then become: 



DDW= ^--\Wf - gcosQZ^2- eF^2 



2 



W, 



D^iy = 0, 
e 



SF12 



d\WY 



-W^^W. 



= cos(9SFi2 - sini95Zi2 + ig' 



sm( 



ddc^) 







4 cos^ [ 



\m\ 



2 cos2 Q 

2\(\)[(\? 



Z(tf 



V 

~2 



2W[d^i- 



9 



► + c.c. . 



(6) 

(7) 

(8) 
(9) 

(10) 

(11) 
(12) 



O. 



d^^A^^ Z^^Wfj^ and their field strengths: O 



2 cos (9 
where 1) = V ^ ig cos 6Z is a covariant derivative. 

In Ref. [M] the equations of motion are solved in the 
Bogomolny limit, Mz = M//, where Mh = V2Xv and 
Mz = gv/{2cos6) are the masses of the Higgs and Z 
bosons, respectively. Here we solve - partially following 
Ref. [H] - the equations of motion for arbitrary mass 
of the Higgs in the region B ^ Be close to the phase 
transition point, \B — Bc\ <C Be. The later condition 
implies that the quantity e = |iy|/MvF <C 1 can serve as 
a small expansion parameter. 

The combination of the two equations of motion, ^ 
and (|8| , with the requirement of the minimization of the 
energy density (Ig]), lead to the simple Abrikosov equation, 

VW^t)W = 0, (13) 




FIG. 2. (a) The superconduct ing W condensate ( |14[ ); (b) the superfluid Z condensate (17); (c) the Higgs expectation value (17); 
and (d) the superconducting (20) and (e) superfluid (21) transport coefficients as the functions of the transverse plane coordi- 
nates xi and X2 at the Higgs mass Mh = l.lMz (the solid surfaces) in the background magnetic field B = l.OlBf^ directed 
along the xs axis. The semitransparent surfaces correspond to Mh = 2Mz in Figs, (a), (b), (c) and (d); and to a^:^ = plane 
in Fig. (d). The red line in Fig. (c) gives the B = Higgs expectation value, = v/V^. Figure (f) illustrates the 3d regions of 
space predominantly occupied by the superconducting electric current J^ of the W bosons and the superfluid neutral flows Js 
of the Z bosons generated by a weak test electric field ^ext parallel to the strong magnetic field B — l.OlBf^ at Mh = l.lMz. 



valid up to corrections of the order of O(e^). This equa- 
tion has nontrivial periodic solutions known as Abrikosov 
lattices. Following Abrikosov [16], we choose a general 
solution of Eq. (IlS]) as a sum over lowest Landau levels: 






\ {\z\^ -\-z^)-7Tu^n^ -\-27Tunz 



(14) 



where Lb = y/27T/{eB) is the magnetic length and u is 
an arbitrary real valued parameter. In order to ensure a 
regular structure of the lattice, the complex coefficients 
Cn are usually chosen in a periodic manner: C^+at = Cn- 

The solution with N = 1 and u = 1 defines the square 
lattice of the original Abrikosov's solution p^. However, 
the global energy minimum is reached for the equilat- 
eral triangular lattice at A/" = 2 with Ci = ±iCo and 
ly = v^/a/2 ^ 0.9306 in agreement with earlier stud- 
ies devoted to the W condensation [14, 15 . The energy 
density ^ is then minimized numerically for fixed val- 
ues of the magnetic field B and the Higgs mass Mh [the 
later is chosen in a phenomenologically interesting range 
of Mh = (1.1 . . . 2.0)Mz]. At B ^ Bf"^ the condensa- 
tion of the W bosons, Fig.jlja), makes the energy density 
smaller compared to its value in the trivial ground state. 
Fig. [lib). The heavier the Higgs boson the weaker the 
effect of the magnetic field on condensates and, as we will 
see, on transport coefficients. 

Equations (fisl) and ^ imply d{B - e\W\'^ /2) = 0, up 



unit cell A of the Abrikosov lattice is 

where the integration constant - given by the integral 
over the unit lattice cell A - guarantees the conserva- 
tion of the magnetic flux, f ,dzdz B{z) = Area(^) • ^ext- 
Thus, the magnetic field ( |15| ) becomes transversally non- 
uniform due to the backreaction of the W condensate. 



Then, the solution of Eqs. (10) and (11), are given by 
the nonlocal expressions for the Z and Higgs condensates: 



Z = Zi^ iZ2 
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to 0(e) terms. The solution carrying a finite energy per 



respectively, with A = Bd = df -\- 82^ The remaining 
equation ([l2| is satisfied automatically up to O(e^). 

In the ground state the W, Z and Higgs condensates 
are functions of the transversal coordinates xi and X2, 
as visualized in Figs. [2|a), (b) and (c), respectively. The 
expectation value of the Higgs field falls down as the 
magnetic field rises, with a slope which becomes weaker 
as the Higgs mass rises. Figs, mc) and^c). 

We point out that the ground state of the vacuum at 
B > Bf^ is a "tandem" phase which is, simultaneously, 
an electromagnetic superconductor and a neutral super- 
fluid. Indeed, introducing an infinitesimally weak test 
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I a superfluid vortex on top of a superconductor vortex] | superfluid vortex | and antivortexj 



FIG. 3. Kaleidoscopic ground state: the density plots of 
the phases (left) of the W condensate (14) and (right) of 
the Z condensate ( |17| ) in the transversal (xi,X2) plane at 
B — 1.015^^. The endpoints of the cuts in the phases 
(shown by the circles) are the superconductor vortices and 
superfluid ( ant i) vortices, respectively. 



electric field E^^^ one can prove - with the use of the 
equations of motion Q - the following transport laws 
for the electromagnetic current J^ = d^Fj^^ ex SC/dA^ 
and the neutral Z-boson current J^ = d^ Zj^^ ex 8C/8Z^ ^ 
respectively: 

d^oJi]ix) = -K^{xi,X2)-Ef\ a[oJ| = 0, (18) 
5[o4(a;) = -K^(a;i,X2)-^r, a[oJf = 0, (19) 

where i = 1,2. The transport parameters for the electro- 
magnetic and neutral currents are, respectively, as fol- 
lows, 



^^{xi,X2)=e^\WWxi,X2), 



K^{xi,X2) 
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(20) 
\W\\xuX2). (21) 



They are functions of the transversal coordinates xi^2, as 
shown in Figsj2|d) and (d), respectively. 

Equation (|18y implies anisotropic superconductivity of 
the ground state at B > Bf^ similarly to an analogous 
phenomenon in QCD [TOl |Tl]: a weak electric field in- 
troduces a resistance-fee growth of electric current which 
continues streaming after the field is switched off. Equa- 
tion (19) implies an anisotropic superfluidity of the neu- 



tral Z currents, and, simultaneously, it illustrates a very 
unusual physical effect: an external electric field induces 
a superfluid current of neutral particles which are flowing 
frictionlessly along the magnetic field axis. 



From the point of view of the electric conductivity 
properties, a ground state of the vacuum can either be 
a superconductor or an insulator due to Lorentz symme- 
try (indeed, a dissipative behavior, like in the Ohm law, 
is inconsistent with the Lorentz symmetry of the vac- 
uum). Thus, the absence of the electric resistance (and 
vanishing shear and bulk viscosities) in the B > Bf^ 
phase are protected by a remnant Lorentz symmetry in 
the (xq^xs) plane. Similar Lorentz-protection arguments 
apply to the superfluid property. 

The superconductivity coefficient ( [2Q| ), averaged over 
the transversal (xi,X2) plane, is a linearly growing func- 
tion of the magnetic field B, Fig. [ijd), at 5 > Bf"^. 
The superfluid coefficient (21) is a sign-changing func- 
tion. Fig. l2Ve), of the transversal coordinates Xi^2 which 
has vanishing mean value if averaged over the transversal 
plane. Thus, an external electric field Ef^^ gives rise to 
(i) a nonzero net electric current along the magnetic field 
axis and (ii) a neutral superfluid space-dependent flow in 
opposite directions with vanishing net current. Fig. p[f). 
The superconductivity and superfluidity always appear 
in a tandem. The transversal components £^^^2* induce 
neither superconducting nor superfluid currents. 

It is known that the ground state of the vacuum at 
B > Bf^ is an equilateral triangular lattice of the vortex 
defects in the W field O |T5] (we call these vortices the 
"superconductor vortices"). At the vortex positions the 
field W ex W^ + iW2 vanishes. Fig. [2|a), and its phase, 
arg(W), winds around the vortex position. We find that 
the ground state has a much more complicated structure 
in the neutral sector: the state has a hexagonal lattice 
of the "superfluid" vortices (and antivortices), character- 
ized by the vanishing Z = Zi -\- iZ2 field. Fig. |2Vb), and 
positive (negative) winding numbers [17 . The combined 
"kaleidoscopic" pattern of the vortex lattices, superim- 
posed on the density plots of the phases of the supercon- 
ducting, W and the superfluid, Z, fields, are shown in 
Fig. [3J All superconductor vortices are always located 
on top of certain superfluid vortices. 

In conclusion, in this paper we have unravelled the 
detailed structure of the electroweak sector of the vac- 
uum in a strong magnetic field and discovered a kaleido- 
scopic pattern of superconducting W^-boson vortices and 
superfluid Z-boson vortices and antivortices. We have 
shown that application of an electric field generates a su- 
perconducting current in tandem with neutral superfluid 
currents. The neutral currents flow predominantly along 
the Z-boson vortices and have vanishing net current. 
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